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We present a theory for Coulomb drag between two mesoscopic systems. Our formalism expresses 
the drag in terms of scattering matrices and wave functions, and its range of validity covers both 
ballistic and disordered systems. The consequences can be worked out either by analytic means, such 
as the random matrix theory, or by numerical simulations. We show that Coulomb drag is sensitive 
to localized states, which usual transport measurements do not probe. For chaotic 2D-systems we 
find a vanishing average drag, with a nonzero variance. Disordered ID-wires show a finite drag, 
with a large variance, giving rise to a possible sign change of the induced current. 
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Moving charges in a conductor exert a Coulomb force 
on the charge-carriers in a nearby conductor, thus induc- 
ing a drag-current (see Fig. Q). This happens whenever 
the distance between the two conductors is of the same 
order as the average distance between charge carriers. In 
recent years Coulomb drag in two-dimensional systems 
has been studied extensively Q and has provided valu- 
able information about the interactions between adjacent 
extended electron gases. 

Coulomb drag of mesoscopic structures has been ad- 
dressed in the case of ID-systems both within the Boltz- 
mann equation approach Q and for Luttinger liquids 
with strong interwire interactions ||. 

The study of fluctuations in the mesoscopic regime was 
recently initiated by Narozhny and Aleiner and it 
was established that fluctuations will dominate at tem- 
peratures smaller than the Thouless energy. This was 
predicted to be the case even for large extended sam- 
ples, such as those used in the 2D experiments [Q. While 
Ref. Q concentrated on structures larger than the phase- 
breaking length, £<p, here we study Coulomb drag of meso- 
scopic samples smaller than 1$. Experimentally there is 
so far only little work on drag in structures with L < 
(<P • We believe this would be an extremely promising 
new direction for the study of mesoscopic transport prop- 
erties, since it gives an opportunity to directly study in- 
teraction and correlation effects in mesoscopic structures. 
Especially disordered mesoscopic systems are known to 
exhibit interesting and unusual physics and the same can 
be expected for disordered Coulomb drag systems - per- 
haps even more so because Coulomb drag in addition to 
the dependence on the transmission properties also has 
a strong dependence on the nature of the wave function 
inside the mesoscopic region. We note that Coulomb cou- 
pling also has interest in other contexts, such as capaci- 
tive coupling of a mesoscopic conductor to the environ- 
ment, charge pumping in quantum dots, or spin-polarized 
transport S. 




FIG. 1. Schematic geometry of a mesoscopic Coulomb drag 
experiment Q. 

Before presenting the technical details we state our 
main results. We develop a formalism for studying drag 
in mesoscopic systems, and apply it to a number of spe- 
cial cases. In the case of ID-wires, studied numerically, 
we find that even a small amount of disorder induces fluc- 
tuations, such that the drag can exceed the ballistic limit, 
be strongly suppressed, or even change sign. The sign 
change is a general feature of mesoscopic drag, which we 
also demonstrate for chaotic systems. Here arguments 
based on random matrix theory show that the drag is 
zero on average, while the fluctuations are finite. The 
zero average drag can thus be taken as a test of the de- 
gree of ergodicity of the system under investigation. Fur- 
thermore, we address the importance of localized states 
in the sample ||. While localized states do not usually 
effect the ordinary transport properties, they turn out 
to be important for the transconductance. The reason 
is that the electron-electron interaction allows for tran- 
sitions in and out of the localized states, which become 
visible at temperatures smaller than the level spacing, 
giving rise to peaks in the transconductance when they 
cross the Fermi level. We also find a temperature depen- 
dence which is very different than the T 2 -dependcncc 
found for extended states. 

General formulation - Using linear response theory 
similar to Refs. P,M we find the Coulomb drag to second 
order in the interaction between mesoscopic subsystems, 
U\2, taking the isolated systems to be otherwise non- 
interacting. The general formula for the dc transcon- 
ductance in the case of two mesoscopic conductors, as 
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illustrated in Fig. 1, is given by 
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where A is the three point correlation function (Ipp), as 
explained in Ref. |jTo|| . Eq. (1) generalizes the results of 
Rcf. J9|,[l0| to systems with broken translation-invariance. 
For the case of mesoscopic conductors it becomes 

Ai(w, r, r') = -2in 2 h ^ 6}{r, v',sp- hu) 



x \n F (ep - huj) - n F (sp)] + (r <-> r'; to — > — w). (2) 
Here 
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where £ Q = e a — e and i labels the subsystem. The ma- 
trix elements are given by = (a\P\j) and p l a g{r) = 
(a\r)(r\P), where |a)'s are the eigenstates of the uncou- 
pled subsystem with energies e a . Using scattering states 
as the basis we get ilg = ■^S eayefj j a/ 3, where the matrix 
j can be expressed in terms of the 2N x 27V scattering ma- 
trix S [0 as j = (t 3 - SMS). Here, r% n , = ±S nn > with 
plus for n belonging to right moving scattering states and 
minus for the left moving states. 

Some general features immediately follow from Eq. (Q). 
The usual cancellation of velocity and density of states, 
which is central in the derivation of the Landauer- 
Biittiker formula, occurs only for /* , whereas for p l a g 
this is not the case. Consequently, in contrast to individ- 
ual subsystem conductances Ga, G 2 \ peaks at the on- 
set of new modes in either of the subsystems. Secondly, 
we notice that the sum over |/3) mixes both propagating 
and evanescent modes. This means that apart from the 
transmission properties also localized states are probed 
by measuring drag conductance. Finally, we notice that 
the outcome of Eq. can have any sign, which is di- 
rectly related to lack of translation-invariance. 

The low temperature limit also follows readily from 
Eq. ([j]). The factor sinh~ 2 cuts off the frequency inte- 
gration and we can expand the A's to lowest order in u>. 
This gives A cx u with the sum over states restricted to 
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those at the Fermi level (fg = Ep — £p): 



A^,r,r') = 4u;TT 2 h 2 lmJ20 P (r,r' ,e F )S(^). (4) 



We immediately sec that the transconductance in this 
limit becomes proportional to T 2 , in accordance with 
the usual Fermi liquid result for electron-electron scat- 
tering. Note however that the low temperature expan- 
sion breaks down when the temperature becomes smaller 



than the level spacing of the discrete, i.e. localized states, 
which we discuss in detail below. At higher tempera- 
tures the T 2 -behavior is replaced by a weaker tempera- 
ture dependence (e.g. for a quasi ID-system, G 2 \ oc T 
for kT > hvp/L as considered in Ref. S). Here we how- 
ever concentrate on the low temperature dependence. 
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FIG. 2. Relative fluctuations of the transconductance as a 
function of the mean free path for the following lengths of 
the two ID-wires: k F L — 100 7r/3 (squares), 200 7r/3 (cir- 
cles), and 300 7r/3 (triangles). The full lines are the results 
of Eq. Q which is shown diagrammatically by the lower left 
inset. The upper right inset shows a typical histogram based 
on 10 4 disorder configurations. 

One- dimensional wires - Next we consider as an il- 
lustrative example two disordered ID-wires, which we 
solve both numerically and analytically using perturba- 
tion theory. The one-dimensional case shows that a small 
amount of disorder can lead to large fluctuations for the 
drag response and even reverse the sign. The reason for 
this is that inter- wire interaction induced forward scatter- 
ing gives rise to a drag response provided it is combined 
with disorder induced backscattering. In contrast, in the 
case of clean wires the backscattering is induced solely 
by the interwire interaction, and therefore the disordered 
case is larger by a factor of order (7l)Ui 2 (0)/Ui 2 (2kF), 

with U 12 (q) = J L dx 1 dx 2 e l ^- x ^U 12 (x 1 ,x 2 ) being 
the Fourier transformed interaction and 1Z is the reflec- 
tion coefficient, which is inversely proportional to the 
mean free path (TZ\ ~ L/£. We can show this explicitly 
by considering the lowest order perturbation theory in 
disorder potential, corresponding to the diagrams shown 
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in Fig. 2, and for long wires k F L S> 1 we find 

{[SG 21 (e)} 2 ) 1/2 [2(^ 1 )(^ 2 )C/ 1 2 2 (2A ;f )c7 1 2 2 (0)] 1/2 
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The denominator is the result G2i(oo) c>c Uf 2 (2k F ) for 
ballistic wires. For the realistic case where U\ 2 (2k F ) <C 
C/i2(0) we see that the fluctuations of the drag can ex- 
ceed the average value. This is in contrast to the fluctu- 

1/2 

ations of the diagonal conductance ([SGu] 2 ) , which 
are vanishing compared to the mean value (Gu) = 
(2e 2 /h)(l - (Hi)) ~ 2e 2 /h in the limit of weak disorder. 
Fig. 2 displays the prediction of Eq. (||) along with the 
numerical results described below and very good agree- 
ment is seen. 

In order to solve the ID model numerically, we study 
Eq. (4) on a lattice using the method of finite differ- 
ences The method offers a way of studying dis- 
ordered systems by ensemble averaging over different 
disorder-configurations ]l2|] . In our numerical example, 
we use a bare long-ranged Coulomb interaction and the 
Anderson model with diagonal disorder |H| . We have nu- 
merically studied the drag as a function of the mean free 
path I and the length L, choosing the Fermi energy cor- 
responding to a quarter-filled band, and for a separation 
given by k F d — 1. We calculate both Gn, G22, and G21. 
Since the potentials in the two wires are uncorrelated 
we in general have Gn 7^ G22, but (Gn) ~ (G 22 ) and 
((5Gn) 2 ) ~ ({5G 22 ) 2 ). Our numerical results for distri- 
butions, mean values, and fluctuations for Gu are in full 



agreement with the results of Abrikosov 14 1. In the de- 
localized regime <»Lwe find as expected that disorder 
has almost no effect on Gu and (Gu) ~ 2e 2 /h with very 

small fluctuations. Fig. || shows ([<5G 2 i(£)] 2 ) 1 ^ 2 normal- 
ized by the drag G2i(oo) in the ballistic regime as a func- 
tion of k F £. The expected 1 /^-dependence is born out by 
the numerical calculations and we also find that the fluc- 
tuations increase with the length of the wires. The inset 
shows a typical histogram of the drag conductance show- 
ing that depending on the disorder configuration G 2 \{€) 
can be either higher or lower than in the ballistic regime. 
Furthermore, note that in agreement with the arguments 
given above the drag conductance shows a sign reversal 
for some disorder realizations. 

Localized states - The low temperature expansion 
Eq. (|4j> , which results in a Independence, is only valid if 
|/3) belongs to a continuum of states. To investigate the 
effects due to localized states we split A in two parts, 



A = Ad + A;, where the first term is given by Eq. (Q) 
while the second term is due to scattering in and out of 
localized states, 

A 1 (x,j,; W ) = -2iA fclWIv) 

/3fElocalizcd 

x [4>(x, y; Bp - huj)(n F (ep - hu>) - n F {epj) 

+4>{y, x; Ef3 + hu)(n F (E/3 + Tluj) - n F (£p))] , (6) 

where the localized states have been chosen to be real 
functions, and 

4>{x, y; e)=^6 (£ q ) S (£ 7 ) ( 7 |I|o) (a\x){y\ 7 ). (7) 

At low temperatures we can approximate 4>{x, y;s ± 
hui) w (f)(x,y;sp), which allows the temperature depen- 
dence to be extracted by integration over ui in Eq. (jl]). 
Furthermore, for temperatures less than the level spac- 
ing the response will be dominated by the coupling to the 
localized level lying closest to the Fermi level. There are 
thus three different types of contributions correspond- 
ing to the response due to localized/delocalized states 

in each subsystem, G 2 i = \ {gd-d + 9l-d + 9l-l) where 
gd-d T 2 . Let us consider, say, gi-d in some detail. We 
find 
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where £1 is the energy of the localized level lying closest 
to the Fermi energy. A similar calculation gives 
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cosh[(ei - e F )/2kT) cosh[(e 2 - £ F )/2kT] 
The relative strengths of these terms can be estimated as 
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with A being the interaction area. Thus at low temper- 
ature the contributions due to localized states will domi- 
nate. The temperature dependence is very different from 
the usual T 2 law, and it may even be temperature in- 
dependent if both £1 and £2 lie on the Fermi level. By 
adjusting the Fermi energy or system parameters, one 
can use the drag response to probe the properties and 
statistics of localized states. 

Random matrix theory - We now discuss the statistical 
properties of the transconductance. This is important in 
order to determine the size of the Coulomb drag for an 
ensemble of disordered mesoscopic systems, such as sug- 
gested in Fig. [I]. Our starting point is the low tempera- 
ture result (|4|) (neglecting localized states) . For the calcu- 
lation we need the statistical properties of the S-matrix, 
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the eigenstates, and the eigenvalues. We assume that the 
region where the subsystems couple by Coulomb interac- 
tions are disordered and that they can be described by 
random matrix theory Jl5| ]. This means that the eigen- 
values and the wave functions are assumed to be uncorre- 
cted and furthermore that the current matrix elements 
I a /3 are uncorrelated with the value of wave functions. 
The latter follows from the fact that the current matrix 
elements are independent of position and may be eval- 
uated outside the disordered region, and hence do not 
correlate with the wave functions inside the disordered 
region. With these approximations 



With typical numbers for GaAs 2DEG structures the 
fluctuations of the transresistance are of the order of 0.1 
Ohm, which should be measurable. 

In conclusion, we have studied drag in mesoscopic sys- 
tems and argue that measurement of transconductance 
provides an interesting new method for investigation of 
the electronic properties of these systems. 

We thank C. W. J. Beenakker and M. Brandbyge for 
useful discussions. 
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The average of the current matrix element is evaluated 
using standard RMT Jlj|, and both with and without 
time reversal symmetry we find (I ai ) = (?i/2to)((t 3 + 
S^r 3 S) ai ) cx r 3 7 , and since the second average in (A) 
is symmetric with respect to interchange of a and 7 
we get (A) = and of course therefore (Cr2i)=0. The 
fluctuations are, however, nonzero and involve the av- 
erage (A(w, r, r')A(o;, s, s')) and hence the combination 
((S^T 3 S) af3 (S^T 3 S) a >p>), which in the limit of a large N 
becomes (2N)~ 2 8 a pid a 'p. Interestingly, again the result 
is not changed by breaking of time reversal symmetry, 
in contrast to the UCF case, where the results with or 
without an applied B- field differ by a factor of 2 |l5| . 
The variance of the A then reads 
(A(u;,r,r')A(£,s,s')) _ C(r, r', s', s) - C(r, r', s, sQ 
tt 2 uju ~ (2iV) 2 

where C is a correlation function involving four density 
matrices 

C(r,r',s,s') = ^2 (p a f3(r)Pp a >{r')p a >0>(s)Pt3>p(s') 



- (2^ (A(r ' 0) (A(S ' S ' }> (A(r ' S ' }> {A{r '" s)) 

to lowest order in 1/kpl. Using the average spec- 
tral function relevant to the 2D case (A(r)) ~ 
(m/2h 2 )exp(—r/2i)J (kFr), and assuming, in addition 
to Uf(. ^> 1 also that I 3> r s , where r s is the screening 
length, we obtain the estimate 
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